† Corresponding author. E-mail:
The thermal properties of pure graphene and graphene–impurity (impurity = Fe, Co, Si, and Ge) sheets have been investigated at various pressures (0–7 GPa) and temperatures (0–900 K). Some basic thermodynamic quantities such as bulk modulus, coefficient of volume thermal expansion, heat capacities at constant pressure and constant volume of these sheets as a function of temperature and pressure are discussed. Furthermore, the effect of the impurity density and tensile strain on the thermodynamic properties of these sheets are investigated. All of these calculations are performed based on the density functional theory and full quasi harmonic approximation.
Graphene is a monatomic layer of carbon atoms in a honeycomb lattice. Figure
The first important conceptual advantage introduced by density functional theory (DFT)[6] is the possibility to describe the ground state properties of a real system in terms of its ground state electronic charge density instead of the far more complicated wave functions. All of the calculations in this paper are performed based on density functional theory using the generalized gradient approximation (GGA) (for the exchange–correlation energy)[7,8] and WIEN2K code.[9] The WIEN2K computational code allows performing electronic structure calculations of solids using density functional theory. It is based on the full-potential linearized augmented plane wave (LAPW) + local orbitals (lo) method, one of the most accurate schemes for band structure calculations. The 4 × 4 × 1 and 2 × 2 × 1 hexagonal honeycomb lattice graphene sheets with 22 Bohr distances between the layers are simulated. By replacing one of the C atoms in the 4 × 4 × 1 and 2 × 2 × 1 graphene sheets with Si, Ge, Fe, and Co impurities, the Si–graphene, Ge–graphene, Fe–graphene, and Co–graphene sheets are created with two different impurity densities (1/32 = 3.1% and 1/8 = 12.5%) (see Fig.
In the harmonic model, the vibrations of atoms in a crystal are treated as 3nN non-interacting phonons with volume independent frequencies ωi, where n and N are the number of atoms per primitive cell and the number of cells in the solid, respectively. The lack of anharmonicity leads to unphysical behaviors,[12] such as zero thermal expansion, infinite thermal conductivity, etc. The harmonic approximation at any crystal geometry (the quasiharmonic approximation (QHA)) is the simplest way of accounting for the anharmonic effects.[13–16] In QHA, the non-equilibrium Helmholtz free energy is
The full quasi harmonic approximation as implemented in the GIBBS code[5] is used to investigate the thermodynamic properties of pure graphene, Si–graphene, Ge–graphene, Fe–graphene, and Co–graphene sheets with different impurity densities.
To investigate the structural properties of pure graphene sheet and graphene sheets with Si, Ge, Fe, and Co impurities, the total energy as a function of volume for these sheets within GGA is calculated and the results are fitted with the Murnaghan equation of states.[17] The results of these calculations for the impurity density of 3.1% are shown in Fig.
The dominant effect on the bulk moduli of these sheets is from the degree of partial ionicity and strong covalency character. The bulk modulus generally increases with increasing covalency. The ionicity effect is to reduce the bonding charge amount and hence to reduce the bulk modulus.
The bulk modulus of the graphene sheet increases when one of its C atoms is replaced by Co, Fe, and Ge impurities. This is due to the increase of hybridization between the 3d orbitals of the Co, Fe, and Ge atoms and the s, p orbitals of their nearest neighbor C atoms. In our previous paper,[18] it has been shown that the impurity-nearest neighbor distance of the Si–graphene sheet is larger than that of graphene sheets with Co, Fe, and Ge impurities. The bulk modulus of the graphene sheet decreases when one of its C atoms is replaced by a Si impurity. This trend is due to the decrease of covalency with increasing nearest neighbor distance and larger Si atomic radius than C atomic radius.
Furthermore, the bulk moduli of these sheets decrease and are in the same order of magnitude when the impurity density increases to 12.5%. We also apply a small strain along the diagonal of the sheet to investigate the thermodynamic quantities at tensile strain. The mechanical properties and shear modulus of graphene have been reported,[19–30] however to the best of our knowledge, there are no experimental or theoretical values of the equilibrium unit cell volume and bulk modulus of this sheet with Si, Ge, Fe, and Co impurities reported.
The thermal properties of pure graphene sheet and graphene sheets with Si, Ge, Fe, and Co impurities are investigated in the temperature range of 0–800 K (below the melting point of pure graphene sheet) and the pressure range of 0–7 GPa. The pressure effect on the volume (EOS curves) at room temperature (299 K), 500 K, and 800 K are calculated and the results are shown in Fig.
A better presentation of these changes is shown in Fig.
The heat capacities of pure graphene, Si–graphene, Ge–graphene, Fe–graphene, and Co–graphene sheets at constant volume (CV) and constant pressure (CP) as a function of pressure are investigated and shown in Fig.
The CV and CP of pure graphene, Si–graphene, Ge–graphene, Fe–graphene, and Co–graphene sheets at tensile strain from 0 to 4 GPa are calculated and compared with the corresponding results under hydrostatic pressure in Fig.
The CV of these sheets with different impurity densities as a function of temperature at different pressures is plotted in Fig.
The Dulong and Petit limit of these sheets decreases by about 1.18% with increasing pressure. The CV of these sheets increases with increasing impurity density at given pressure and temperature. The CV of these sheets decreases slightly with increasing pressure. The coefficient of volume thermal expansion (α) for pure graphene and graphene sheets with Si, Ge, Fe, Co impurities as a function of temperature at different pressure is calculated and shown in Fig.
The volume of pure graphene sheet and graphene sheets with Si, Ge, Fe, and Co impurities increases with increasing temperature at a given pressure. The bulk modulus of pure graphene sheet at T = 0 K increases when one of the C atoms in the graphene sheet is replaced by Ge, Fe, and Co impurities and decreases when it is replaced by a Si impurity. The bulk modulus of these sheets decreases with impurity density increasing. The bulk modulus of these sheets under tensile strain increases compare to the corresponding value under hydrostatic pressure.
With increasing pressure and temperature, the bulk modulus of these sheets decreases and tends to a constant value. CV is proportional to T3 at low temperature and tends to the Dulong and Petit limit at high temperature. CV and CP of the strained sheets significantly decrease compare to those of these sheets under hydrostatic pressure at low temperature. α increases rapidly at low temperature (between 100 K and about 200 K) and tends to increase linearly at high temperatures for a given pressure.
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 | |
11 | |
12 | |
13 | |
14 | |
15 | |
16 | |
17 | |
18 | |
19 | |
20 | |
21 | |
22 | |
23 | |
24 | |
25 | |
26 | |
27 | |
28 | |
29 | |
30 |